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Abstract. - We investigate the effect of dipolar interactions in one-dimensional systems in connec- 
tion with the possibility of observing exotic many-body effects with trapped atomic and molecular 
dipolar gases. By combining analytical and numerical methods, we show how the competition be- 
tween short- and long-range interactions gives rise to frustrating effects which lead to the stabiliza- 
tion of spontaneously dimerized phases characterized by a bond-ordering. This genuine quantum 
order is sharply distinguished from Mott and spin-density wave phases, and can be unambiguously 
probed by measuring non local order parameters via in-situ imaging techniques. 
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Introduction. — Cold atom gases confined in reduced 
4iniensionality represent an ideal system to observe many- 
body phenomena related to the prominent role played 
Ijy quantum fluctuations [HH]. Recent experimental ad- 
vances have paved the way to the investigation of quantum 
magnetism, notable examples being the demonstration of 
super-exchange interactions in bosonic gases [3] , the time- 
evolution of spin impurities [11[5] , the observation of frus- 
Ijrated classical dynamics [S] and the engineering of Ising 
Hamiltonians [7]. 

' New opportunities in this direction are now stimulated 
by the prominent progresses in cooling and controlling 
ultracold gases of magnetic atoms and polar molecules, 
which provide tunable platforms where the effect of long- 
fanged dipolar interactions can play a dominant role in 
determining the many-body dynamics [51-[T^. Such pro- 
gresses have opened a new door for the investigation of lat- 
tice models beyond the conventional Hubbard paradigm, 
where long-range interactions can compete with local ones 
on the way to unveil richer many-body physics |13H15j . 
Much attention has been devoted up to now to extended 
Bose-Hubbard models, where a new phase of matter, the 
Haldane insulator (close analog of the Haldane phase in 
spin-1 chains) has been predicted to occur [TMTO] . How- 
ever, not much is known on other possible magnetic phases 



for fcrmionic dipolar gases in optical lattices [15], which 
present close analogy to the so called extended Hubbard 
model (EHM) [MSO]- 

In this Letter we show that one dimensional (ID) Hub- 
bard models with long-range interactions support a non- 
trivial insulating phase characterized by bond-order [231 
US [571 ISO] [a bond-order density- wave (BOW)] due to 
the competition between dipolar and on-site interactions. 
This phase has attracted notable interest in recent years 
in the context of strongly correlated electron systems. Its 
existence is now well-established but it has been longly 
debated in a series of theoretical studies [^MHO] . The op- 
portunity of realizing and observing such states of matter 
may shed new light on a series of issues, from its dynamical 
properties to its ground state robustness. In order to prove 
its existence in dipolar Hubbard models, we combined an- 
alytical and state-of-the-art numerical methods based on 
the density-matrix-renormalization-group (DMRG) algo- 
rithm |31I32) . Remarkably, the effect of dipolar interaction 
is evident already in the weak-coupling approach, leading 
to the prediction of a BOW phase within one-loop order. 
This is in sharp contrast with standard EHMs - on-site and 
nearest neighbor interaction only-, where a simple formal- 
ism is instead unable to capture its existence [121 ■ While 
the detection of the BOW phase is in general challenging 
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Fig. 1: Panel a): numerical phase diagram for dipolar mix- 
tures described by the Hubbard Hamiltonian in Eq. ^ and 
corresponding cartoon picture for the different phases. A 
region with dominant bond-order develops close to the line 
U ~ V 2 ■ Errors in determining phase transition points 
are smaller than symbol sizes. Panel b): qualitative behav- 
ior of different string correlation functions Oc,3 in the different 
phases (CDW, SDW and BOW, from top to bottom): in the 
BOW phase, both order parameters display long-range-order. 



due to the limited extension in parameter space and its 
lack of density-like order parameters, we show how the re- 
cently developed in situ imaging techniques [33134] provide 
an ideal route toward the unbiased identification of such 
phases since the BOW phase is uniquely identified by the 
long-range order of non-local parity correlations [351136] . 

Spin-1/2 dipolar Fermi gas. Dipolar particles 

confined in a one dimensional tube can be described by 
the following microscopic Hamiltonian: 



H = f dxY, |v4(^) 

•^ CT=1,2 '- 



2m 



-dl 



ila(x) 



^ / dx dy\p„(x){y{x,y) ^U5{x,y))pa'{y)\, 



(1) 



cr.cr' 

where m is the particle mass, V'ct(V'o-) are cre- 
ation(anmhilation) operators of cr ~ 1,2 fcrmions, and 
Pa = 'iply'i'a- The short-distance contribution of inter- 
species interactions lA can be tuned by using Fcshbach res- 
onances or confinement induced, while the dipolar interac- 
tion V is controlled by means of electric (polar molecules) 
or magnetic (magnetic atoms) fields [TS]. Once confined 
on a sufficiently deep optical lattice, an effective descrip- 
tion in terms of Hubbard Hamiltonian leads to the discrete 
formulation (see, e.g., [T]): 

H = -y^t(c|^^Cg-,i+i +h.c.) + ?7y^ni^,ri.2,i 

a.i i 

+ ^^y<y,<T'i\i- i\)na,tna',j, (2) 



where the on-site interaction U is now given by a com- 
bination of local and dipolar potential, and the last term 
describes long-range dipolar interactions, Va^a'{\i ^ i\) = 
V / {\i — i\f' , which we consider purely repulsive from here 
on. In the equal mass, equal interaction case, the sys- 
tem inherits a global SU{2) symmetry, which is preserved 
by the long-range tail, and is reduced to U{1) for gen- 
eral parameter choices. From here on, we will focus on 
the former case and consider a balanced half-filled chain, 
Ni = N2 = N/2 = L/2, where N^ is the number of par- 
ticles in the spin state a and L the system size. Like in 
the case of the EHM, the phase diagram of Eq. ^ is de- 
termined by the competition between local and non-local 
interactions. Before searching for specific BOW instabili- 
ties, we illustrate the competing mechanism in the atomic 
limit i = 0, as such competition will then lead to spin- 
frustration at the origin of bond-order itself. 

For dominant U interactions, the system ground state 
is a spin density wave (SDW) . On the other hand, a dom- 
inant dipolar interaction will minimize the energy by im- 
posing double occupancies every second site, thus stabiliz- 
ing a fully gapped charge density wave (CDW) . These two 
ground states, illustrated in the insets of Fig. |TJ become 
energetically degenerate along the line Uc = 3C(3)F/2 (^ 
is the Riemann zeta function), which determines a phase 
transition between SDW and CDW. Deeply in the quan- 
tum regime t ~ y ~ [/, however, quantum fiuctuations 
may enhance the emergent frustration close to the clas- 
sical transition line and lead to a different critical sce- 
nario. For the EHM, where only nearest-neighbor interac- 
tions are considered, it was argued by Nakamura [33] that 
an additional phase with dominant charge bond-order in- 
stability occurs between the SDW and the CDW phases. 
This phase is characterized by a spontaneous spin-Peierls 
dimerization, manifest in a charge polarization on alter- 
nating bonds and by the formation of spin dimers on the 
bonds, and constitutes a notable example of dimerization 
in strongly correlated systems due to frustration. A 2D 
analog of this phase has been recently discussed for spin- 
less dipolar fcrmions in layers [37] . 

The existence and extent of the BOW phase in the 
EHM have been intensively debated mainly because i) 
the instability is not captured by one-loop g-ology calcu- 
lations based on bosonization (while more refined meth- 
ods recently provided analytical evidence of it f^l^ ). 
and ii) numerical results were not consistent due to the 
small extent in parameter space of such a phase, and to 
the difficulty of providing accurate location of the critical 
line [331|5S1[171|3D]. In the following, combining analyti- 
cal and numerical methods, we show that dipolar systems 
support spontaneous spin-Peierls dimerization, and that 
the corresponding BOW phase occupies a larger region in 
parameter space with respect to the EHM. First of all, we 
show that one-loop g-ology is sufficient to establish the ex- 
istence of bond-order in the weak-coupling regime. Then, 
the existence and extent of the BOW phase are bench- 
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marked with DMRG simulations. 

Low-energy field theory. — We now present a qual- 
itative study of Eq. ^ within the bosonization frame- 
work [38j . As a first step, we express the fermionic lattice 
operators in terms of continumir chiral fields ipr,/T{x): 



R/Ly 



J 



V^K.(^)e*'"'+V'L(a^)e- 



ikj^xi 



(3) 



with X = ja, a the lattice spacing and kp ~ T:N/{2La) the 
Fermi momentum. We then apply the standard bosoniza- 
tion mapping introducing density and phase fiuctuation 
fields ip^, {}„ for the two species in order to map the orig- 
inal fermionic problem onto a bosonic one: 



'0fl,<T 



VR.ct ,;[^^-^^1 



V'l 



VL, 



jW„+-d„] 



(4) 



where 77^,0- are the so-called Klein factors. We proceed by 
first bosonizing independently the single species Hamilto- 
nians, and then consider the coupling between them. This 
allows us to partially keep track of the long-ranged dipo- 
lar interactions by using known results on single channel 
Luttinger liquids (LL) with dipolar interactions [39ti42] . 
The single species Hamiltonian can then be bosonized as 
follows: 



H„ 



hv 
Y 



dx 



K 



Kid^-d,) 



g / (ia;cos[v47r(^(j] 



(5) 



where v is the sound velocity, and the Luttinger parameter 
K is weU approximated a.s K = {1 + 1.46n * V/t)-^/"^ [42] . 
An additional mass term due to the lattice is also present. 
While these results strictly hold in the continuum, they 
well capture the quantitative dependence of K with re- 
spect to interaction and filling fraction even in lattice mod- 
els |43] . We can now proceed and treat local and non-local 
interspecies interactions as a perturbation on the top of 
Hi + H2. After introducing density and spin collective 
fields: 

ipi ± ip2 „ ^1 ± l?2 



Vc/s 



1} 



(6) 



^/2 ' ''' V2 ' 

the Hamiltonian can be decoupled into a density and spin 
sector, H ^ He + Hs, where the physics in each sector is 
described by a sine-Gordon model: 



He 



hv, 



c/ s 



dx 



Kr/s 



+ X,/,(9,i9,/,)2 



+ 9c/s I dx cos[y ^iTi If c/ 



(7) 



Notice that terms coupling density and spin degrees of 
freedom are also present; in the following we will neglect 
their effects, which are supposed to be small away from 
the regime V,U <^ t due to their large scaling dimen- 
sion. Within first order g-ology, the parameters of the 
effective Hamiltonian can be estimated as follows: first, 



we compute the effect of the interspecies interactions on 
the quadratic part of the Hamiltonian by approximating 
the interspecies dipolar interactions up to next-nearest- 
neighbor. The effective Luttinger parameters then read: 



K. 



c/ S 



IK 



K 



± 



U + W/A 

47rf 



(8) 



The coefficient gs of the mass term in the spin sector is 
then fixed by imposing SU(2) spin symmetry in Hg. This 
leads to a Berezinskij-Kosterlitz-Thouless (BKT) transi- 
tion [35] which takes place at Kg ~ 1, that is, in the 
weak-coupling limit, Vc ^ ~ C//2.34. In the repulsive 
regime [/, 1/ > 0, the density sector is instead expected 
to be always gapped since Kc < 1, except on the line 
gc — which identifies a Gaussian transition (at strong 
coupling; this line can be unstable toward Akp mass terms 
when Kc < 1/4). Since the coefficient of the mass term is 
gc oc {U + 7V/4:), the Gaussian transition can be identified 



at Vc' 



(G) 



C//L75. 



The phase diagram can be subsequently mapped out 
by considering the dominant orders as in the case of the 
EHM [23j. FotV < C//2.34, the system is in a SDW phase, 
with a gapless spin sector and dominant correlations of the 
form {{ui^^ ~ ni4)(ni+j;,-f- - n^+x^C))- Oi^ t^e other hand, 
a (mass) density-wave is formed above the critical value 
V > Vc , where the dominant correlations are of the form 
{{jii^^ + ni,x){ni+x,-\ ^ rii+x.x)) ■ In the intermediate regime 
y{BKT) ^y ^ y{G) ^ neither SDW nor CDW order are 
stable, and the system exhibits a BOW, characterized by 
both a finite spin and density gap, and a dominant order 
described by the parameter 



m-^^Y. 



_a-Ci+l,CT + Cj_|_]^^g.Ci_CT 



(9) 



This indicates that dimers are spontaneously formed on 
nearest-neighbor bonds, a phenomenon usually called as 
spontaneous spin-Peierls dimcrization. 

As such, we conclude that, even within one-loop g-ology, 
the low-energy theory already predicts a finite region of 
parameter space where bond order is stable. This is due 
to the effects of the long-range dipolar tails, which consis- 
tently affect the Luttinger parameter of the single species 
Hamiltonians, and lead to different conditions for the BKT 
and Gaussian transition. However, the treatment is sup- 
posed to work only in the weak-coupling regime, and the 
provided estimate on the Gaussian line does not indeed 
capture the entire dipolar interactions. In order to con- 
firm this qualitative view, we provide in the following a 
numerical analysis of Eq. ^ . 

Numerical results. — The accurate determination of 
the phase diagram of the system is a very challenging task. 
As far as the simpler EHM is concerned, despite the great 
effort put in its numerical study over the last decade, few 
results have shown good agreement with each other [Mll^ 
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[571[3D| . Mindful of the difficulties in the determination of 
the phase diagram for such a system, we carefully checked 
our methods and codes on the EHM with both periodic 
(PBC) and open boundary conditions (OBC); discarded 
weights of the order 10 ~^ allow us to estimate the phase 
transitions within 1% with the most recent and accurate 
results [3D]. 

We calculate several physical quantities to determine as 
accurately as possible the two phase boundaries of Eq. ([2|), 
keeping in the DMRG simulation terms up to |j — i| = 5 in 
the off-site dipolar interaction part. We verified for small 
system sizes that the inclusion of longer range terms does 
not affect significantly the quantities of interest. 

CDW-BOW transition. Starting in the CDW phase, 
as the dipolar interaction strength is decreased the system 
enters in a BOW phase. As a first probe of the transition, 
we calculate the density gap 

Ac= lim [E{N + 2,0) + E{N-2,0)^2E{N,0)]/2 

(10) 
where E{N, Sz) is the ground state energy oi N — L par- 
ticles with total magnetization Sz — (N^- — Ni)/2. Due 
to the competition between long range and on-site interac- 
tions the gap has to take a minimum value at the transition 
point. To locate the phase boundary as precisely as pos- 
sible we also calculate the BOW order parameter defined 
as 

(B) = UmL^^\{BL/2 - Bl/2+i)1 (11) 

where Bi is the operator in Eq. ^. {B) is the amplitude 
of the oscillation of the BOW operator at the center of the 
chain, defined in such a way that Friedel oscillations are 
weaker [44]: a non- vanishing value of {B) will be a clear 
signature of the BOW phase. 

A further signature of the phase transition is given by 
the LL parameter Kc- For a periodic chain it can be ex- 
tracted from the (density) static structure factor 



Sc{q) 



l^e»?('.^-0((„^„,) -{nk){ni)) 



(12) 



k,l 



with q = 27r/L. Within LL theory K^ = limq_>.o 7:Sc{q)/q 
is finite only on the phase transition line for a continuous 
transition (while is always zero instead for first order tran- 
sitions). Since we are dealing with a finite size system we 
expect to see a sharp peak at the transition line varying 
the dipolar interaction. 

In order to get the correct values in the thermodynamic 
limit for the quantities described above, a careful finite 
size scaling analysis must be carried out. The density gap 
Ac is extrapolated fitting the data with a fourth-order 
polynomial in 1/L. For the BOW order parameter (B) 
we assume that in the centre of the chain the amplitude 
of Friedel oscillations is proportional to l/L^" [33] to ex- 
trapolate its thermodynamic value. 

We run several DMRG simulations for different dipolar 
interaction strengths V and system sizes (till L = 128) 
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Fig. 2: Panels a) and b): Luttinger parameters K^ and Ks 
respectively, as functions of V for different L; orange triangles 
are the values of Ks in the thermodynamic limit (TDL). Panel 
c): BOW order parameter (B) for various V. Straight lines 
are guides for the eye, whereas dashed lines fit the numerical 
data using the appropriate scaling laws. 



with OBC. We keep up to r?i = 1256 states and per- 
form 6 sweeps to calculate the density gap and the BOW 
order parameter. In this way the truncation error, i.e., 
the weight discarded in the renormalization procedure, 
is at most of order 10~^, while the error on the ener- 
gies is approximately 10~®. To calculate the static struc- 
ture factor (|12p . we instead use PBC on smaller systems 
{L = 16, 24, 32, 40, 48), but the same number of states and 
sweeps as for OBC. 

The results for the CDW-BOW transition are reported 
- ior U = 4t - in panels a) and c) of Fig. [5] and panel a) 
of Fig. [3J In particular one sees that Kc develops a peak 
(sharper as the system size increases) when V/t ~ 2.4, the 
BOW order parameter (B) vanishes as V/t — 2.44 while 
it is still finite for V/t ~ 2.4 and the density gap has a 
minimum also at V/t = 2.4. Thus we can infer that the 
transition between the CDW and the BOW phases occurs 
for V/t ~ 2.4. 

SDW-BOW transition. As the strength of the dipo- 
lar interaction is further decreased, Peierls dimerization is 
destroyed and the system enters in a SDW with a uniform 
distribution of the density and no gap in the spin sector. 
The BKT nature of this transition makes its location chal- 
lenging when evaluating the spin gap 



A, = lim [E{N, 1) + E{N, 0)] 



(13) 



since it is exponentially small close to the transition line 
|38j . A valid alternative is provided by the spin-static 
structure factor 



Ss{q) 



k.l 



Mk-l) 



Mst 
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(14) 
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Fig. 3: Panels a) and 6j: charge Ac and spin gaps As respec- 
tively, for various V; a magnification is illustrated in the inset 
in panel b). Panel c): density and spin (x25) gaps in the ther- 
modynamic limit for various V. Straight lines are guides for 
the eye, whereas dashed lines fit the numerical data using the 
appropriate scaling laws. 



Indeed LL theory predicts that Kg = liiriq^Q 'KSs{q)/q is 
zero in the spin gapped phase and Kg = 1 in the gap- 
less one. Logarithmic corrections prevent Kg to reach the 
latter value even for long chains. Nonetheless, such correc- 
tions have been shown to vanish in the frustrated J — J' 
model when the system forms dimers [45] . The same is 
true for the EHM [37], since the BOW-SDW transition 
should have the same nature, and thus also for our model. 
In the BOW phase Kg = near the transition only for 
very large systems. Following [57] we estimate the transi- 
tion point when, at fixed U, Kg crosses 1 from above as V 
is increased. The thermodynamic limit of the spin gap A^ 
is obtained, as for the density gap, fitting the data with 
a fourth-order polynomial in 1/L, reproducing the spinon 
band structure near the band edges. 

The results of the DMRG simulations arc reported in 
panel h) of Fig. [2] and panel h) of Fig. |3] First of all one 
can see that Kg crosses 1 for V/t < 2.1 and an accurate 
finite size scaling allows us to locate the transition point 
for V/t = 2.089. This result for the phase transition point 
is also confirmed by both {B) and the spin gap, that have 
a small, but still finite value (precisely {B) = 5.1 x 10^'^ 
and As = 1.6 x lO^'^) at V/t = 2.08. 



Dipolar 2-leg ladder. — Wc have also analyzed 
a different configuration, which requires single species 
fermionic dipoles. The gas is confined in two parallel tubes 
at a distance g and the dipoles (the external electric or 
magnetic field) form an angle 9 with the plane of the wires 
as sketched in Fig. 2] a). The system can be modelled by 




Fig. 4: Panel a): schematic cartoon of ladder configuration. 
Panels b)-c): Kc and Ks as a function of the angle 9 for differ- 
ent system sizes for the case di = d2 — 1.5 and g = 1.0 



the Hamiltonian 



H, = - 



E ^ [4. 



E E^^ 

a,fe— 1,2 i.j 



Ci+l,a 



1 — 3 sin 



E 'i 



,2 ^i,a^j,a 



\^-3\ 



~'^i,a^j,b 



(15) 



where a,b = 1,2 are the tube indices, t the hopping pa- 
rameter and d the dipole momentum, and r denotes the 
distance between particles in different tubes. We per- 
formed DMRG simulations with PEC for system sizes 
L ~ 12, 20, 24, 32 to calculate the LL parameters Kc and 
Kg for several values of the angle 6. The results are re- 
ported in panels b) and c) of FigdJ by applying an anal- 
ysis similar to the one presented above, we conclude that, 
whilst a BOW phase may indeed exist for intermediate 
tilting angles, its extent, indicated by the maxima of the 
Luttinger parameters, would be much smaller than that 
for the single tube Fermi mixture. 

Experimental regimes and probes. Polar 

molecules and magnetic atoms offer strong versatility in 
tuning interactions. In the latter case, the ratio V/U can 
be independently tuned by means of Feshbach resonances, 
which have been already reported for bosonic ^^Cr [Tl] and 
^^^Er [12] isotopes. In case of molecular gases stable un- 
der two-body recombination, accurate estimates of on-site 
interactions will be required. An alternative approach can 
employ Feshbach molecules of strongly magnetic atoms 
such as Er or Dy, effectively increasing by a factor ^ 8 the 
strength of the dipole-dipole interactions [12] • 

The existence of the BOW phase can be indirectly 
probed spectroscopically as follows: first, the density gap 
is estimated by means of lattice modulation spectroscopy, 
indicating the onset of the Gaussian transition. Subse- 
quently, the spin gap is estimated by means of RF spec- 
troscopy. The BOW phase can be located in the intermedi- 
ate region between the two transition point, as illustrated 
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in Fig. [3J:. A more solid, direct probe of the exotic na- 
ture of the insulating state is the long-range nature of the 
parity order parameters: 



e+x 



i+x 



Osix) = {Y[: 



MS' 



Ocix) = (H ' 



j=e 



j=i 



[3] 
[4] 
[5] 
(16) [6] 

[7] 
[8] 



which are well-defined order parameters for general Hub- 
bard models [3S] . As discussed in Ref. P^[55] and as illus- 
trated in Fig.\^b)-d), such correlation functions properly 
distinguish the CDW, SDW and BOW phases. In partic- 
ular only the latter phase has long-range order in both Os 
and Oc- Although measuring parity order parameters is 
a difficult task, recently Oc has been measured for cold 
atoms with short range interaction loaded in an optical 
lattice by in situ imaging of the many-body wave function 
via atom fluorescence [23]. 

Conclusions. — In summary, we have provided a de- 
tailed study of how dipolar fcrmionic mixtures support 
exotic insulating states with dominant bond-order. We 
have underpinned the corresponding Hubbard phase dia- 
gram by combining numerical and analytical techniques. 
Remarkably, differently from the standard EHM, we find 
that the BOW phase can be found already at the level of 
one-loop g-ology. This is due to the long-range nature of 
the interaction which enhances the frustration of the sys- 
tem. We also briefly discuss how single species loaded in 
ladder geometries do not represent a favorable platform to 
investigate bond-ordered phases. 

Experimentally the non-trivial correlations embodied in 
the bond-density- wave phase can be faithfully captured by 
spin and density string order parameter by means of in 
situ imaging. This makes cold atoms in optical lattices 
an ideal setup for the investigation and demonstration of 
bond-order instabilities in strongly correlated systems. 
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